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Abstract 

We exploit the quantum coherence between pair-produced and in ip{3770) decays to make 
a first determination of the relative strong phase differences between Kg-K~^-K~ and D^' 

Kgir'^TT^ , which are of great importance in determining the CKM angle 7/03 in B~ — ^ D^{D^)K~ 
decays. Using 818 pb~^ of e^e' collision data collected with the CLEO-c detector at = 3.77 
GeV, we employ a binned Dalitz-plot analysis of X^vr'^^vr" and K^t:~^t:~ decays recoiling against 
flavor-tagged, CP-tagged and K^'k^tx~ tagged events to determine these strong phase differences. 



* Deceased 
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I. INTRODUCTION 



A central goal of flavor physics is the determination of all elements of the CKM matrix [1] , 
magnitudes and phases. Of the three angles of the b — d CKM triangle, denoted a, /3, and 
7 by some, 02, <Pi, and 03 by others, the least- well determined is 7/^3, the phase of 
relative to Vcb- It is of great interest to determine 7/^3 using the decay — > K^D^, 
since in this mode, the 7/^3 value obtained is expected to be insensitive to new physics 
effects in B decays. Here, is either D'^ or D'^, and both decay to the same final state, 
and so their amplitudes add. Sensitivity to the angle 7/^3 comes from the interference 
between two Cabibbo-suppressed diagrams: b — > cus, giving rise to B- K'D^, and the 
color and CKM suppressed process b ucs, giving rise to B- K-D^. One of the most 
promising D'^ decays for measuring 7/^3 using this method is D° — > i^^TT+TT", because it 
is Cabibbo favored (CF) for both and decays, thus providing large event yields. To 
make use of this decay, however, the interference effects between B~ K~D^[-^ K^t:^t:~) 
and B^ K^D^{-^ K'gH^n) need to be understood. These interference effects can be 
understood and measured using CLEO-c data. 

We first write the amplitude for the B^ decay as follows: 

A{B^ ^ K^b\b^ ^ X°7r+7r-(x,y)) oc /^(x, y) + rBe^^±/s(x, y). (1) 

Here, x = fn^o^+i y = "^^-o^- are the Dalitz-plot variables in the decay, 

fD{x,y){fQ{x,y)) is the amplitude for D^{D^) decay to Kg7i~^7i~ at {x,y), tb is the ra- 
tio of the suppressed to favored amplitudes, and 9± = 7- where 5b is the strong phase 
shift between the color-favored and color-suppressed amplitudes. Ignoring the second-order 
effects of charm mixing and CP violation [2, 3], we have fr){x,y) = fD{y,x), and Eq. 1 can 
then be rewritten as: 

A{B^ ^ K^b',D' ^ Kln^n-ix,y)) oc Mx,y) + tbc''^ foiy, x). (2) 

The square of the amplitude clearly depends on the phase difference ASd = SD{x,y) — 
6D{y,x), where 6D{x,y) is the phase of fo^x^y). Thus, for the determination of 7/^3, one 
must know Ado^x^y). 

Previous analyses extracted A.5r){x,y) by fitting a fiavor-tagged — > XgTr+TT" Dalitz 
plot to a model for decay involving various 2-body intermediate states [4-6]. Such an 
approach introduces a 7° ~ 9° model uncertainty in the value of 7/^3, ^ which would be a 
limiting uncertainty for LHCb [7] and future i?-factory experiments. 

In the analysis presented here, we employ a model-independent approach to obtain 
A5d(x, as suggested by Giri et al. [2], by exploiting the quantum coherence of — 
pairs at the ■0(3770). Because of this quantum correlation, K^-k^-k^ and K^-k^ti^ decays 
recoiling against flavor tags, CP-tags, and — > K^ti^ti~ tags, taken together provide di- 
rect sensitivity to the quantities cosA^^i and sinA^^. This measurement will result in a 
substantial reduction in the systematic uncertainty associated with the interference effects 
between B' K-D'^{-^ X^+tt") and B' K-D%-^ X^tt+tt). 



^ BaBar claims 5° uncertainty on 7/^3 in [5] by combining KgK~^K and Kg7r~^TT modes. 
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II. FORMALISM 



Giri et al. proposed [2] a model-independent procedure for obtaining A5d{x, y), as follows. 
The Dalitz plot is divided into 2jV bins, symmetrically about the hne x = y. The bins are 
indexed from —i to i, excluding zero. The coordinate exchange x ^ y thus corresponds 
to the exchange of bins i ^ —i. The number of events in the i-th bin of a flavor-tagged 
X^TT+TT" Dalitz plot from a decay is then expressed as: 

Ki^ADj^fDix,y)\^dxdy^ADFi, (3) 

where Ad is a normalization factor. The interference between the and amplitudes is 
parameterized by two quantities 

Ifoix, y) Wfoiy, x) I cos[A5d{x, y)]dxdy, (4) 



^/FW- 
and 

Si = ^^= f Ifoix, y)\\fD{y,x)\sm[A6D{x,y)]dxdy, (5) 

where the integral is performed over a single bin. The parameters q and Sj are the amplitude- 
weighted averages of cos A^^ and sin A5/j over each Dalitz- plot bin. It is important to note 
that Ci and Si depend only on the D° decay, not the B decay, and therefore these quantities 
can be measured using CLEO-c data. In principle they could be left as free parameters in a 
Kg7r~^7r~ Dalitz-plot analysis from decays, but their values can be more precisely 
determined from correlated D^D^ pairs produced in CLEO-c. 

Though the original idea of Giri et al. was to divide the Dalitz plot into square bins 
[2], Bondar et al. noted [8] that increased sensitivity is obtained if the bins are chosen to 
minimize the variation in A^^ over each bin. Thus, wc divide the Dalitz phase space into 
J\f bins of equal size with respect to ASd as predicted by the BaBar isobar model [4]. In the 
half of the Dalitz plot m'^{KgTT~^) < w?{Kg'n~), the i^^ bin is defined by the condition 

27r(i - Z/2)/N < A5d{x, y) < 27r(i - l/2)/N, (6) 

The — i*'* bin is defined symmetrically in the lower portion of the Dalitz plot. Such a binning 

with A/" = 8 is shown in Fig. 1. One might suspect that because we arc using a model to 
determine our bins, we are not free of model dependence. In fact any binning is correct 
in that it will give a correct, unbiased answer for 7/^3, at the cost of larger uncertainties 
compared to an optimal binning with respect to A^^. 

We now describe how CLEO-c data can be used to determine q and Sj. The event yields 
in the i^^ bin of both flavor-tagged and CP-tagged Z^*^ K^-k^ti^ Dalitz plot arc required. 
Because the ■;/'(3770) has C = — l,_the CP of the Kgii^Ti^ decay can be determined 

by reconstructing the companion in a CP eigenstate. With a CP-tagged D° — > Kgn~^TT~ 
decay, the amplitude is given by: 

fcp±ix,y) = -^[fDix,y)± fDiy,x)], (7) 
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for CP-even and CP-odd states of a — > Kg-K^n^ decay. Since the event rate is propor- 
tional to the square of this amphtude, the number of events in the i^^ bin of a CP-tagged 
Dahtz plot is then: 



Mt = hcP±{K, ± 2c,^K,K., + K.,), (8) 

where hcp± = /2Sf is a normalization factor that depends on the number, S"/, of single 
flavor-tagged signal decays, and the number, S'^, of single CP-tagged signal decays. Thus, 
access to q is enabled by measuring the number of events, Mf, in a CP-tagged i^'^vr+TT" 
Dalitz plot, and the number of events, Ki, in a flavor-tagged K^Tt^-K^ Dalitz plot. 

Unfortunately, as evident from Eq. 4, the sign of A^^i is undetermined in each of the i 
bins. However, sensitivity to both q and Si can be obtained by analyzing — > K^tt^tt~ 
vs. K^n^Ti' data. The amplitude for ?/'(3770) decays to two i^'^vr+vr" decays is as 

follows: 

f{x,y,x,y) = -j= . (9) 

The prirned and unprimed Dalitz-plot coordinates correspond to the Dalitz-plot variables of 
the two D'^ K'^TT+TT- decays. Defining Mij as the event rate in the i^^ bin of the first and 
the j*'^ bin of the second —>■ K^ti^ti^ Dalitz plots, respectively, we have: 



Mij = Korr{KiK^j + K^.Kj -2,/K~K~K~K'j{ciCj + SiSj)). (10) 

Here, hcorr = ^ds/2>S'|, where Nj;)j^ is the number of DD pairs, and as before Sf is the 
number of flavor-tagged signal decays. Equation 10 then relates the product (qCj -|- SiSj) to 
the measured yields of events in the flavor-tagged — > i^'^vr+vr" Dalitz plot [Ki^s) and the 
yields in the K'^tt+ TT VS. — > Kg7i~^7T (Mjj's) Dalitz plots. The sensitivity to this 

product leads to a four- fold ambiguity: change of sign of all Cj or all Sj. In combination with 
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the CP-tagged analysis though, where the sign of q is determined, this reduces to a two-fold 
ambiguity. One of the two solutions can be chosen based on a weak model assumption [4]. 

The decay — > K1ti~^ti~, due to its close relationship with — > KgTT^7T~ , can be used 
to further improve the Cj and Sj determination. Since the Kg and mesons are of opposite 
CP, and we assume the convention that A{D^ A'^tt+tt") = A{D^ — >■ Kg7r~7r^), it then 
follows that A{D^ KItt+tj-) = -A{D^ Kln-n+). Then, for irOTr+vr-.^the Dalitz-plot 
rates of Eq. 8 (CP vs. /sr^+TT") and Eq. 10 (D° ^ i^^+TT" vs. i^^TT+Tr") 

become: 

Mt = hcP±{K[ T 2c',^KiK'_, + K'_,), (11) 
Mij = h^r[KiK'_j + K_iK'^ + 2^KiK'_jK_iK'j{cic'j + s^s'^)], (12) 

for CP vs. D° — > X^TT+TT" and K^gT^^i:' vs. Ki^7r~^7r~, respectively, where c^, 

are associated with D° — > X^tt+tt" decay. 

For — > ii'°7r"'"7r~ decays to benefit our determination of q and Sj, we must determine 
the differences Acj = c' — Cj, and Asj = s'- — Si. In addition to the relative sign change in 
Eq. 11 and Eq. 12, doubly Cabibbo suppressed decays (DCSD) of D^/D^ also contribute 
with opposite signs in — > Kgn'^n' and D° — > K^Tr'^n' decays. We can see this by 
inspecting the decay amplitude for each Dahtz plot 

A{Kln+n-) = -^[A{K\+n-) + A{K\+Tr-)], (13) 
v2 

AiKlir+TT-) - ^[A(X°7r+7r-) - A(XV7r-)]. (14) 
v2 

The effect of this relative minus sign is to introduce a 180° phase difference for all DCSD K* 
resonances in the K^ii^n^ model. We can use U-spin symmetry to relate the amplitudes 
for resonances of definite CP eigenvalue, e.g. Kg j^p^{770). We find that these states acquire 
a factor of re"^. To convert a D° ^ X°7r+7r~ model to the corresponding X^tt+tt" 
model, we multiply all DCSD amplitudes by —1 and multiply each CP cigenstate amplitude 
by (1 — 2re*'^), with r = tan^ 9c and 5 = 0° fixed for every resonance (here 6c is the Cabibbo 
angle). We then determine central values for the corrections, Aq = c- — q, Asj = s- — Sj using 
the Kg7r^7r~ BaBar model [4]. We ascribe uncertainty to both the choice of r and 

5, as well as the usage of the BaBar model to determine the uncertainties on Aci and Asj. 
The former are estimated by varying the phase S between and 27r, and r by ±50%. For 
the latter, we estimate Aq and Asj using the Belle [6] and CLEO [9] Kg7r~^7i~ isobar 

model fits, and take the largest resulting deviation from the value found with the BaBar 
model as a model-dependent systematic uncertainty. The total systematic uncertainties in 
the corrections are the quadrature sum of these two uncertainties. The central values and 
uncertainties on Aq and Asj are shown in Table I. 

III. EVENT SELECTION 

We analyze 818 pb~^ of e+e~ collision data produced by the Cornell Electron Storage Ring 
(CESR) at E^^ = 3.77 GeV and collected with the CLEO-c detector. The CLEO-c detector 
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TABLE I: Predicted values for Acj and Asj with the systematic uncertainties. 



Ac,; 



0.099 ± 0.040 -0.034 ± 0.068 
0.167 ± 0.029 -0.064 ± 0.084 
0.327 ± 0.122 -0.013 ± 0.097 
0.253 ±0.192 0.133 ±0.136 
0.077 ± 0.061 0.041 ± 0.080 
0.220 ± 0.084 -0.038 ± 0.065 
0.416 ±0.160 0.095 ±0.063 
0.184 ± 0.024 0.015 ± 0.086 



is a general purpose solenoidal detector which includes a tracking system for measuring 
momentum and specific ionization {dE/dx) of charged particles, a Ring Imaging Cherenkov 
detector (RICH) to aid in particle identification, and a Csl calorimeter for detection of 
electromagnetic showers. The CLEO-c detector is described in detail elsewhere [10]. 

Standard CLEO-c selection criteria for tt^, K^, 7r°, and Kg candidates are used, and 
are described in Ref. [11]. To distinguish electrons from hadrons, we use a multivariate 
discriminant [12] that combines information from the ratio of the energy deposited in the 
calorimeter to the measured track momentum (E/p), ionization energy loss in the tracking 
chamber (dE/dx), and the ring- imaging Cherenkov counter (RICH). For Kg decays, we 
select candidates with |M(7r+7r" -) - Mkq\ < 7.5 McV/c^, and require the decay vertex 
to be separated from the interaction region with a significance greater than two standard 
deviations (except for — > KgTT~^7T~ vs. 5° KgTT~^TT~ candidates). Reconstruction of 
77 — > 77 proceeds analogously to 7r° 77, with the requirement that |M(77) — < 42 
MeV/c^. We form u vt+tt^tt^ candidates and require their mass to be within 20 MeV of 
the nomimal uj mass [13]. 

In this analysis, we reconstruct mesons in several fiavor-tagged modes, CP-tagged 
modes, and in Kgn'^n". From these selected events, we also reconstruct the companion 
from the ^jJ{3770) decay in either Kgn'^n^ or _ft'°7r"'"7r~ to form "double-tags". The single 
tags yields enter our analysis through the Sf and S"^ factors, whereas the double-tags provide 
the Ki, Mi and Mij yields across their respective Dalitz plots. The double-tagged events we 
consider are shown in Table II (all the notations include charge conjugate if not otherwise 
specified.). We thus consider flavor tags: K~7i~^, K'tt'^tt'^, K~t:'^t:'^t:~; semileptonic tag: 
K-e+v- CP-even tags: K+K-, tt+tt", K^t^^tt^, K^tt^; and CP-odd tags: K^7r°, K^gV, 
with tt"/?] — 77, u — >■ 7r"'"7r~7r'', and Kg tt+tt". We also reconstruct double-tag (DT) 
events with D° — >■ Kg ^tt"'"tt~ vs TT+TT as discussed in the preceding section. We do 

not reconstruct tt+tt" in some DT modes when there are two missing particles {K~e'^i', 
and Ki^n^ cases) or the backgrounds are large (as for X^tt^tt", and Kgu). 

The event yield in the i*^ Dalitz-plot bin of each tagged —> i^'^j^ tt+tt" sample is 
determined by evaluating the phase difference for each data point according to the BaBar 
isobar model. The contribution of each isobar to the total amplitude is evaluated as a 
function of all three invariant mass-squared combinations computed directly from the four- 
momentum of the D daughters as described in Ref. [14]. The phase difference is well defined 
beyond and continues smoothly across the kinematically allowed Dalitz-plot boundary as 
shown in Fig. 2. A small number of candidate events (~l-3% depending on tag and signal 
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mode) included in this analysis are reconstructed outside the kinematically allowed region 
due to finite detector resolution. 



TABLE II: Reconstructed Double Tag modes. 



Mode 


Kgvr+vr- Kl 


TT+vr 


Flavor Tags 


K-TI+ 


X 


X 




X 


X 




"TT^ X 


X 


K-e+v 


X 




CP-Even Tags 


K+K- 


X 


X 


7r"'"7r~ 


X 
X 
X 


X 


CP-Odd Tags 




X 


X 




X 

X 


X 




X 


X 



8 



A. Single Tags 



The V^(3770) resonance is below threshold for DDn production, and so the events of 
interest, e^e~~ — > ■0(3770) — > DD, have D mesons with energy equal to the beam energy and 
a unique momentum. Thus, for identifying candidates, we follow Mark III [15] and define 

two kinematic variables: the beam-constrained candidate mass, Mbc = \JeI/c^ — Po/c^, 
where is the candidate momentum and Eq is the beam energy, and Ai? = — i?o, 
where Ej:, is the sum of the candidate daughter energies. Candidate tags are required to 
have A.E within about 3 standard deviations of zero [16]. 

For events with a K~t^^, K^K~, and tt+tt" single-tag (ST) that have no additional 
charged particles, we apply additional selection requirements to suppress cosmic ray muons 
and Bhabha events. We do not allow tracks identified as electrons or muons to be used in the 
tag. We demand evidence of the other D by requiring at least one electromagnetic shower 
in the calorimeter above 50 MeV not associated with the tracks of the tag, where a single 
minimum ionizing particle deposits the equivalent of 200 MeV. For K^K~ ST candidates, 
additional geometric requirements are needed to remove doubly radiative Bhabha events 
followed by pair conversion of a radiated photon. We accept only one candidate per mode 
per event; when multiple candidates are present, we choose the one with smallest |Aii^|. 

The resulting Mbc distributions are shown in Fig. 3. Each distribution is fit to a signal 
shape derived from simulated signal events and to a background ARGUS [17] threshold 
function. The ST yield is given by the area in the signal peak in the mass region from 
1.86 < Mbc < 1-87 GeV. 



B. Double tags with Kgir'^ir 

We form Kgn'^n" DTs by combining a K^n'^n' tag with a ST candidate. We choose 
one DT candidate per mode per event with M closest to the measured D° mass, where 
M =[M{D^) + M{D^)]/2. 

Since the —>■ KgTT~^TT~ vs. — > Kgn~^7i~ sample plays a key role in extracting Sj 
values, we drop the requirement on the flight distance signiflcance for Kg candidates to 
mcreasc the statistics. We find 421 D° K'^tt+ 71 VS. D° Kgir'^ir candidates which 
include about 9% background. We increase the yield by about 15% (additional 54 candidates, 
~15% background) by reconstructing the i^^Tr+vr" vs. /T^tt+tt" candidates when one tt^ 
is not reconstructed. The presence of the tt"^ is inferred from the missing four-momentum 
calculated from the well known initial state and the reconstructed particles. 



C. Double tags with K e^v 

Candidate K'e^u vs. A'^tt+tt^ DTs are reconstructed by combining a K^tt^tt~ ST 
candidate with a kaon candidate and an electron candidate from the remainder of the event. 
Events with more than two additional tracks (aside from the K^t:'^t:~ daughters) are vetoed. 
Signal discrimination for D — > K'e^v uses the variable U = Emiss — c|p„j,s.s|- where Emiss 
and Pmiss are the missing energy and momentum in the semileptonic meson decay, 
calculated using the difference of the four-momenta of the tag and that of the K~ and e"*" 
candidates. For correctly identifled events, U — 0, since only the neutrino is undetected. 
After all selection criteria are applied, multiple candidates are rare for K~e'^u. The U 
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FIG. 3: Data M^c distribution for various tag modes. The solid lines show the total fits, and the 
dashed lines show the background shapes. 



distribution for D K'e^u candidates is shown in Fig. 4. The points with error bars 
are data and the shaded histogram represents a simulation of the background, which is less 
than 1% in the signal region, \U\ < 50 MeV. 



D. Double tags with K£7r° 

The K^T[^Ti~ vs. K^T[^ DT mode is reconstructed with a missing mass technique since the 
K^j^ mesons produced at CLEO-c are not reconstructed. A fully reconstructed Kgn^-n' ST 
is combined with a 7r° candidate, and we compute the recoil-mass squared against the ST-7r° 
system, M^^^^. Signal K^ti^ decays are identified by a peak in M^^^^ at M|,o. Backgrounds 

from D — > -ft'^TT*^, 7r°7r°, and r^vr^ are suppressed by vetoing events with additional unassigned 
charged particles, or ^ 77 or 7r° 77 candidates. We further suppress backgrounds by 
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FIG. A: U = Emi ss — c|pmiss| distribution for K e^v in events with a — > K^tt^tt signal can- 
didate. The points are data and the shaded histogram represents a simulation of the backgrounds. 
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FIG. 5: M^^gg distributions for K^tt^ when one D is identified as Kg-K^n^. Shaded histogram 
represents simulation of the backgrounds. The enhancement of data relative to simulation ~ 
0.9 GeVVc"^ corresponds to the decay i^*°7r° K£7r°7r° where both the and vr*^ from the K*^ 
are undetected. The Dalitz-plot model of this process is not implemented in our simulation. 



making requirements on the energy of showers in the calorimeter that are not associated 
with the decay products of the Kg7i~^n~ or the 7r°. We compute the angle, 9, between each 
unassigned shower and the direction of the missing momentum. For cos 6' < 0.9, we require 
the energy of showers, Eghower < 100 MeV for any single shower. If 0.9 < cos^^ < 0.98, 
we require Eshower < 100 + 250 x (cos 6* - 0.9) MeV. The M^^^^ distribution for K^vr^ is 
shown in Fig. 5. The points with error bars are data and the shaded histogram represents 
a simulation of the backgrounds. Signal candidates are required to be within the range 
0.1 < M^i,, < 0.5 GeVVc^ 



E. Double tags with K^7v~^7v 

Candidate i^^vr+vr" decays are reconstructed in DTs using a similar missing mass tech- 
nique as described in Section III D. We require the signal side (associated with the /^^Tr+Tr" 
candidate) to have exactly two charged tracks. Backgrounds are reduced by applying 7r°, 
1], and Kg vetoes. Using the measured momenta of the tagged and the two additional 
pions, we compute the missing momentum and missing energy on the signal side. We apply 
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FIG. 6: M^j^^ distributions for K^tt+vt" for flavor tags (left), CP-even tags (middle), and CP-odd 
tags (right). The points with error bars are data and the shaded histograms represent simulations 
of the backgrounds. 



the same requirements to the energy of the unassigned showers as described in Section HID. 
The M^j^gg distributions for i^^^yr+TT" are shown in Fig. 6. The points with uncertainties are 
data and the shaded histograms show a simulation of the backgrounds. Signal events are 
required to have a missing-mass squared in the region 0.21 < M^^^^ < 0.29 GeV'^/c^. 



F. Yields in Data 

The ST yields for the tag modes and DT yields for Kg^^n~^7i~ versus different tags 
are shown in Table III. To determine the — > K~e~^i' and —>■ K^ir^ ST yields we 
use the integrated luminosity, measured D^D^ cross-sections [11] and measured branching 
fractions [18, 19]. Combining all modes of the same CP, we show in Fig. 7 the Dalitz-plot 
distribution of CP-even and CP-odd tagged_i5'^ Kgn^-K^ decays. Figure 8 shows the 
corresponding distributions for CP-tagged K1ti^7i~ decays. The clear absence of 

a p^Kg component (CP-odd) in CP-odd tagged Kgir'^n' decays is an illustration of the 
quantum correlations that exist in the '?/'(3770) — > D^D^ decay. For i^T^Tr+yr", p°i^° is 
absent in the CP-even tagged samples. 

The signal-to-background ratios in our ii'^^^Tr+vr" DT samples range from 10 to better 
than 100, depending on tag mode. The tag side Ai?, Kg and u sidebands are used for 
combinatorial and non-resonant background subtraction. On the signal side, the background 
level is 1.9% for K^tt^tt~ after applying the flight significance requirement. This part 
of the background is considered as a systematic error. The background-to-signal ratio for 
the K^vr+vr" signal side is about 5%, of which about 2% is a peaking background from 
i^gTr+vr", — >■ 7r°7r° decays that pass the K^ir^n' selection criteria. We estimate this 
peaking background yield using K^tt+tt" data and a misidentification rate determined from 
a quantum-correlated Monte Carlo simulation. The combinatorial background contribution 
is estimated using the M^j^^ sidebands. The expected yields from these two background 
sources are subtracted from the observed signal yields to obtain background-corrected yields. 

For the Kg7r^7r~ vs. Kgir^n' sample, no Kg flight significance requirement was applied, 
resulting in a_ background-to-signal ratio of ~9%. About 7% (of 9%) of this background 
comes from 7r+7r^7r+7r^ faking Kgn^Ti^. This background is subtracted using a 

Monte Carlo simulation of this decay, where the 7r+7r^7r+7r~ Dalitz-plot structure is taken 
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TABLE III: Single tag and K"/r7r+7r" double tag yields. 



Mode 


ST Yield Kgvr 


"•"vr yield K^n 


"•"TT yield 




Flavor Tai 


gs 




K-TT+ 


144563 ± 403 


1447 


2858 


K-tt+tt' 


258938 ± 581 


2776 


5130 




220831 ± 541 


2250 


4110 


K-e+i^ 


123412 ± 4591 


1356 




CP-Even Tags 


K+K- 


12867 ± 126 


124 


345 


7r'*"7r~ 


5950 ±112 


62 


172 




6562 ± 131 


56 




Kin' 


27955 ± 2013 


229 




CP-Odd Tags 


K'sTt' 


19059 ± 150 


189 


281 


2793 ± 69 


39 


41 


K'su 


8512 ± 107 


83 




K's-K+TT- 




475 


867 



from the FOCUS experiment [20]. The impact of the remaining ~1.9% of background on 
our nominal fit results is small and included in the systematic uncertainties. 

The reconstruction efficiency is defined as the ratio of reconstructed events to generated 
events in each bin. The reconstruction efficiencies are calculated from large Monte Carlo 
samples generated according to the amplitude description of Eqs. 7 and 9 for different 
tag modes. Dividing the observed yields in each 5d bin by this efficiency, we obtain the 
efficiency-corrected yields, Mf" and M^j. 

IV. EXTRACTION OF a AND Si 

We determine the coefficients q, Si by minimizing the negative log- likelihood function 

-21og£ = -2Y.\ogP{Mt,{Mt))^cp,Ki.^.-) 

i 

-2^1ogP(M±, {Mt)\cp,Kl.+.-) 

i 

-2^1ogP(M,„(M,,)) 

-2 ^l0gP(M,,-, {Mij))^K^^,+,-,Kln+.-) 

+x^ (15) 

where (Mj^) is calculated according Eqs. 8 and 11, and (Mj^) is calculated according Eqs. 10 
and 12, and P(M, (M)) is the Poisson probability to get M events given the expected 
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FIG. 7: CP-even tagged i^'^vr+vr Dalitz plot (a), and its ^.^(Tr+Tr ) projection (b). CP-odd tagged 
KgTr~^7r~ Dalitz plot (c), and its m^(7r+7r~) projection (d). 



number, (M). In our nominal fit, a penalty term 



Ci - Aci 



Asi 



»\2 



SAsi 



(16) 



constrains c[ and s'^ to difi^er from Cj and Si. respectively, by their expected difi^erences Ac,, 
Asi, within errors. (Those errors, 6Aci, 6Asi, are the systematic uncertainties shown in 
Table I.) This constraint has little impact on q but is important for Si and will be relaxed 
and tightened as a systematic variation. 

From Monte Carlo studies, we found that DCSD decays in flavor tag modes {K^n'^, 
K^n^n^, K^n^n^n^) lead to a signiflcant bias in the K^'^^s due to an interference of the 
wrong flavor of the — * K^^j^ti^ti^ decay; this results in a signiflcant bias in the values 



of cl'^ and s\'' from the — >■ i^^7r'''7r vs. — > Kg^^Tr^Ti analyses. Therefore, for the 
— > K^n'^n' vs. D° — > X^Tr+Tr" analysis, we use only the D° — > K~e'^u tagged K^tt+tt" 



sample and the ^ K tt"*" tagged K^tt'^t: sample for counting S^^ ^ and K\ ' yields.^ For 



All three hadronic flavor tag modes are used to determine k\'^ for the CP-tag vs. -ftT^ ,^7r+7r determina- 



S/L' 
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FIG. 8: CP-cven tagged i^'^vr+vr Dalitz plot (a), and its rn^(7r"'"7r ) projection (b). CP-odd tagged 
i^£7r"'"7r~ Dalitz plot (c), and its m^(7r+7r~) projection (d). 

the latter, we estimate the biases and adjust the values using the correction factor: 



\Ado 



Hcrc r = \A{D^ — * K~^n )/A{D^ — > K 7r+)| and Skw are the ratio of amplitudes of the 
DCSD to CF decay and the relative strong phase, respectively. The amplitude ratio squared, 

= (3.44 ± 0.01 ± 0.09) x 10"^ and Sxn = (22 ± 16.3)° are taken from Ref. [16]. This 
correction factor is estimated in each of our eight Dalitz-plot bins using the BaBar D° — > 
KgTT^7T~ Dalitz-plot fit amplitude [4]. The model dependence of this correction is negligible. 
Uncertainties on these corrections due to the uncertainty on Sktt are small and are included 
in our systematic uncertainties. 

The fitting procedure was tested using a simulated C-odd D^D^ Monte Carlo sample 
where we performed 100 toy K^n^n^ vs. Kg7i~^7r~ experiments with Cj and Sj taken from 
the BaBar model. The means and widths of the pull distributions of the q and Sj parameters 



tion of c. 



(') 
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TABLE IV: Fit results for Sj (with Aq 
constrained). See Table I for Acj and Asj. 



and Asi fixed), q, c[, Si and s'^ (with Acj and Asj 



Acj, Asj fixed 



Aci,Asj constrained 



0.742 ± 0.037 

1 0.606 ±0.071 

2 -0.008 ± 0.063 

3 -0.529 ±0.101 

4 -0.889 ± 0.049 

5 -0.742 ± 0.066 

6 0.108 ±0.074 

7 0.403 ± 0.046 



0.004 ±0.160 
0.014 ±0.215 
0.581 ±0.190 
0.138 ±0.217 
-0.053 ± 0.183 
-0.317 ±0.187 
-0.836 ±0.185 
-0.410 ±0.158 



0.743 ± 
0.611 ± 
0.059 ± 
-0.495 ± 
-0.911 ± 
-0.736 ± 
0.157 ± 
0.403 ± 



0.041 0, 
0.072 0, 
0.077 0. 
0.114 -0. 
0.053 -0. 
0.070 -0, 
0.092 0, 
0.046 0. 



840 ± 0.041 
779 ± 0.072 
250 ± 0.078 
349 ± 0.135 
793 ± 0.057 
546 ± 0.080 
475 ± 0.094 
591 ± 0.048 



0.014 ±0.166 
0.014 ±0.216 
0.609 ±0.188 
0.151 ±0.225 
-0.050 ± 0.189 
-0.340 ±0.194 
-0.827 ±0.190 
-0.409 ±0.158 



-0.021 ±0.164 
-0.069 ±0.219 
0.587 ±0.188 
0.275 ± 0.232 
-0.016 ±0.192 
-0.388 ± 0.200 
-0.725 ±0.196 
-0.374 ±0.169 



were consistent with zero and one, respectively, indicating no bias and proper estimation of 
statistical uncertainties. 

To enable the separation of the statistical uncertainty on Cj and Si from the systematic 
uncertainty on Aq and Asj we perform a likelihood fit to (cj, Sj) with the values of (c^, s'j) 
fixed according to Table I. The results of this fit and of the nominal likelihood fit to (cj, Si), 
(c^, s'i) are shown in Table IV. The (statistical) correlation matrix among Cj and Sj in the 
constrained fit is shown in Table V. Note that in Table IV, wc choose the Sj and s- signs 
based on BaBar isobar model predictions [4] to resolve the two-fold ambiguity discussed in 
Section II. 



TABLE V: Correlation Matrix for the q and Sj parameters. Labels 1-8 represents Ci — Cs and 9-16 

reprcsoul .si — .s>^. 



i 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 15 16 


1 


1.000 




























2 


-0.028 


1.000 


























3 


-0.007 


-0.011 


1.000 
























4 


0.035 


0.011 


0.002 


1.000 






















5 
6 


0.073 
0.016 


-0.022 
0.069 


0.006 
-0.003 


0.004 
0.003 


1.000 
-0.104 


1.000 


















7 


0.018 


0.013 


-0.020 


0.005 


0.033 


0.016 


1.000 
















8 


-0.020 


-0.028 


0.008 


0.005 


0.050 


0.013 


0.015 


1.000 














9 


0.024 


0.006 


-0.072 


-0.006 


0.014 


0.014 


0.113 


0.040 


1.000 












10 


0.000 


-0.033 


0.017 


-0.007 


-0.003 


0.038 


-0.001 


-0.001 


-0.060 


1.000 










11 


0.006 


0.007 


-0.025 


-0.002 


0.003 


-0.008 


0.041 


0.029 


0.323 


-0.154 


1.000 








12 


0.004 


0.005 


-0.020 


0.035 


0.002 


0.000 


0.014 


0.005 


0.149 


-0.124 


0.244 


1.000 






13 


0.001 


-0.001 


-0.003 


-0.014 


-0.008 


-0.072 


-0.007 


-0.004 


0.158 


-0.107 


0.340 


0.070 


1.000 




14 


-0.011 


-0.014 


0.078 


0.005 


-0.004 


-0.042 


-0.086 


-0.021 


-0.448 


0.085 


-0.213 


-0.124 


-0.275 


1.000 


15 


0.009 


0.008 


-0.053 


-0.003 


0.004 


0.002 


0.061 


0.013 


0.373 


-0.139 


0.314 


0.228 


0.269 


-0.405 1.000 


16 


0.004 


0.004 


-0.056 


-0.004 


0.003 


-0.007 


0.026 


0.041 


0.234 


-0.096 


0.521 


0.176 


0.243 


-0.133 0.106 1.000 



V. SYSTEMATIC UNCERTAINTIES 



Systematic uncertainties on (q, Sj) and (c^, s'j) come from many sources. Table VI and 
Table VII summarize the main contributions of the systematic uncertainties for Cj and Sj, 
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respectively. Table VIII and Table IX summarize the main contributions of the systematic 
uncertainties for c[ and s^, respectively. 

TABLE VI: Systematic uncertainties for q. 



Cl 


C2 


C3 


C4 


C5 


C6 


C7 


C8 



^ statistics error 0.010 0.015 0.016 0.019 0.009 0.015 0.018 0.008 

Momentum resolution 0.008 0.015 0.012 0.019 0.011 0.012 0.017 0.009 

Efficiency variation 0.004 0.007 0.011 0.008 0.005 0.008 0.010 0.006 

Single Tag yields 0.006 0.007 0.013 0.011 0.005 0.008 0.015 0.008 

Tag side background 0.007 0.007 0.014 0.013 0.006 0.008 0.014 0.011 

fsT^TT+TT- background 0.001 0.002 0.009 0.027 0.012 0.006 0.003 0.002 

Kl-K+T^- background 0.006 0.018 0.004 0.024 0.017 0.012 0.020 0.006 
Multi-Candidate selection 0.002 0.003 0.003 0.008 0.004 0.006 0.003 0.002 

Non-D/£) 0.010 0.016 0.004 0.007 0.004 0.005 0.006 0.005 

DCSD 0.009 0.012 0.005 0.013 0.014 0.010 0.013 0.006 

Sum 0.022 0.037 0.031 0.052 0.032 0.030 0.042 0.021 

TABLE VII: Systematic uncertainties for Si. 

Sl S2 S3 S4 S5 Se S7 Ss 

kP statistics error 0.031 0.027 0.039 0.030 0.023 0.023 0.033 0.026 

Momentum resolution 0.018 0.035 0.023 0.033 0.023 0.022 0.022 0.018 

Efficiency variation 0.018 0.012 0.019 0.010 0.013 0.018 0.013 0.012 

Single Tag yields 0.005 0.001 0.005 0.004 0.003 0.003 0.005 0.003 

Tag side background 0.004 0.001 0.001 0.003 0.001 0.004 0.002 0.001 

K'^TT+TT- background 0.005 0.008 0.030 0.023 0.005 0.003 0.016 0.014 

KItt+t:- background 0.050 0.022 0.018 0.035 0.006 0.024 0.005 0.025 
Multi-Candidate selection 0.036 0.018 0.033 0.012 0.022 0.026 0.028 0.011 

Non-D/D 0.005 0.003 0.004 0.005 0.005 0.005 0.003 0.002 

DCSD 0.023 0.004 0.030 0.019 0.015 0.006 0.027 0.020 

Sum 0.077 0.055 0.076 0.069 0.045 0.052 0.060 0.050 



In the global fit, the fitter does not take the statistical uncertainties associated with flavor 
tagged samples into account. We estimate this part of the uncertainties by varying the input 
variables {K^ ) one by one according to their statistical uncertainties, and by making new 
fits. At the end, we take the quadratic sum of all the variations as the systematic error. 

Since our Dalitz-plot binning results in bins with unusual shapes and in some cases 
very narrow regions (see Fig. 1), the migration of events from one bin to another bin may 
bias our result. The position of an event in the Dalitz plot depends on its momentum 
determination. The systematic error associated with momentum resolution is studied by 
smearing the momentum of a fully simulated Monte Carlo 200 times, according to the 
CLEO detector momentum resolution. The distributions of the results for (q, Si) and {c[, 
s'^) are then fitted with Gaussian functions, and the widths of the distributions are taken as 
the systematic uncertainties. 

The systematic uncertainties associated with Kg^j^ir'^Tr" finding cancel under the assump- 
tion that the efficiency systematic uncertainties are uniform across the Dalitz plot. Under 
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TABLE VIII: Systematic uncertainties for c^. 



Cl C2 C3 C4 C5 C6 C7 Cs 

iff' statistics error 0.009 0.014 0.011 0.013 0.009 0.012 0.017 0.008 

Momentum resolution 0.008 0.016 0.014 0.020 0.010 0.013 0.019 0.010 

Efficiency variation 0.006 0.006 0.009 0.009 0.004 0.007 0.011 0.006 

Single Tag yields 0.005 0.007 0.007 0.007 0.003 0.005 0.006 0.007 

Tag side background 0.005 0.007 0.007 0.007 0.004 0.005 0.005 0.009 

K%'K+-K- background 0.002 0.001 0.005 0.012 0.007 0.003 0.001 0.002 

KI-K+-K- background 0.008 0.020 0.017 0.046 0.019 0.014 0.015 0.006 
Multi-Candidate selection 0.003 0.003 0.001 0.004 0.003 0.004 0.002 0.002 

Non-£)/£» 0.012 0.018 0.005 0.006 0.004 0.004 0.008 0.005 

DCSD 0.011 0.013 0.007 0.013 0.014 0.013 0.014 0.006 

Sum 0.023 0.039 0.029 0.057 0.029 0.028 0.036 0.021 

TABLE IX: Systematic uncertainties for s[. 

Cl C2 C3 C4 Cs C6 C^ C% 

kP statistics error 0.030 0.028 0.037 0.026 0.022 0.025 0.034 0.028 

Momentum resolution 0.022 0.042 0.031 0.042 0.029 0.028 0.031 0.033 

Efficiency variation 0.018 0.012 0.015 0.009 0.012 0.018 0.013 0.013 

Single Tag yields 0.005 0.001 0.005 0.003 0.004 0.003 0.004 0.003 

Tag side background 0.003 0.000 0.002 0.002 0.002 0.004 0.003 0.001 

K^TT+TT- background 0.006 0.009 0.028 0.021 0.004 0.003 0.015 0.015 

Kln+n- background 0.052 0.022 0.013 0.033 0.003 0.025 0.004 0.025 
Multi-Candidate selection 0.038 0.018 0.031 0.007 0.020 0.027 0.028 0.013 

SAci, SAsi 0.034 0.048 0.001 0.006 0.009 0.072 0.050 0.048 

Non-D/D 0.005 0.003 0.004 0.005 0.004 0.005 0.003 0.002 

DCSD 0.021 0.005 0.028 0.016 0.013 0.007 0.028 0.022 

Sum 0.080 0.060 0.072 0.067 0.046 0.056 0.065 0.059 



this condition Eq. 8, Eq. 10, Eq. 11, and Eq. 12 have the same dependence on efficiency. To 
account for a small non- uniformity, we generate a large number of toy experiments where we 
randomly distribute the efficiency of each bin according to a Gaussian distribution (width 
is taken as 0.02) and repeat this process for many times. The widths of the resulting distri- 
butions for (cj. Si) and (c^, s^) are taken as systematic uncertainties. 

The systematic uncertainties associated with the single tag yields are evaulated by re- 
peating the fit with the input values varied by their own uncertainties. Assuming the 
contributions are uncorrelated, we sum in quadrature to obtain the uncertainties on (q, Si) 
and (c^, s'j) due to single tag yields given in Table VI - Table IX. 

The systematic uncertainties due to the estimation of the tag side background are studied 
mode by mode, and the quadratic sum is given in Table VI - Table IX. 

Though we used AE and M(7r+7r~) mass sidebands for the tag side backgrounds sub- 
traction, we did not apply any background subtraction for Kg7r~^7r~ signal side, which is 
believed to be small since we require the decay vertex of Kg to be separated from the in- 
teraction region with a significance greater than two standard deviations. The background 
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level in the signal region is estimated from M(7r+7r^) sidebands. We found there is about a 
1.9% background in the signal region. The systematic uncertainties due to this part of the 
background are estimated using Quantum Correlated Monte Carlo samples. We estimate 
the background contributions from Quantum Correlated Monte Carlo samples, then make a 
new fit with the background subtracted. The differences of the results between the nominal 
fit and the new fit are taken as the systematic uncertainties. 

The systematic uncertainties due to the /T^tt+tt" background shape are considered by re- 
peating the fit assuming the background across the Dalitz plot is uniform. The uncertainties 
due to the estimation of the background level are neghgible. The systematic uncertainties 
due to fiavor-tagged, CP-even tagged, CP-odd tagged, and KgTr~^7i~ tagged Kln~^7r~ sam- 
ples are considered separately and summed in quadrature in Table VI - Table IX. 

It is possible to select a wrong combination when there are multiple signal candidates 
in an event, especially for Kgir'^Tr" vs. Kg^j^n~^Tr~ samples, since there are many pions 
with similar momenta. The systematic uncertainties are studied by applying correction 
matrices to the yield matrices Mij. The corrections are typically 2% (5%) for the Kg7r~^'K~ 
vs. KgT:~^7T~ {KgT:~^7T~ vs. K1tt~^7t~) event samples. 

Monte Carlo simulated continuum events are checked for non.-D°/D^ backgrounds. No 
significant peaking background is seen for double tagged Kg7r'^Tr~ samples. The contribu- 
tions for K~^K~, 7r+7r~, and Kgf] tagged K^-k^ti' samples are also negligible. For other 
samples, there are 1~2% contributions depending on the tag mode. A systematic study is 
performed by assuming the background is uniformly distributed over the Dalitz plot. 

For the DCSD effect, we made corrections to K~t^^ vs. K^ j^t^^t^' yields in Section IV 
by using results from Ref. [16]. The systematic uncertainties due to r is negligible since it 
is precisely measured. The systematic uncertainties due to the strong phase 5 are studied 
by varying it according to its error. For K~t{^t:^ , and i^~7r"'"7r"'"7r~ tag modes, there are no 
relative strong phase measurements, so we consider four cases, 5 = (0°, 90°, 180°, 270°), and 
take the maximum variations among the four cases as the systematic uncertainty. 

The total systematic uncertainties on (q, Sj) and (c^, s'^) ~ excluding the systematic 
uncertainty on Acj and Asj that relate the K^-r^-r^ and K^-k^ti^ Dalitz-plot models — are 
obtained from the quadrature sum of these systematic uncertainties, are shown in Tables VI 
- Table IX. 

In the global fit, Aq and Asj are constrained using a term. The errors on Aq and 
Asj are determined by comparing BaBar, Belle, and CLEO II — > K^ti^ti~ Dalitz-plot 
fit results. The constraint on c and c' can be removed with little impact on the result. 
The constraint on s and s' can be relaxed to a factor of 4, but cannot be removed entirely, 
otherwise, the fit does not converge. To assess our sensitivity to this constraint we consider 
the following 1) we relax the constraint by a factor of 2, i.e. increase the errors by a factor of 
2 and re-fit the data. 2) we fix Aq and Asj and re- fit the data (see Table IV). The maximum 
difference for each (cj, Sj) and (c^, between these fits and the nominal fit is interpreted 
as the systematic uncertainty. An alternate assessment of this systematic uncertainty is the 
difference in quadrature of the errors reported for the "fixed" and "constrained" fits reported 
in Table IV. The quadrature average of these two methods is reported as the third error on 
(ci,Si) and (c^,s9 in Table X and Table XI, respectively. 
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FIG. 9: Comparison of Aq = c • — q between CP-tagged i^^^j^vr+vr CLEO-c data (circles) and 
predictions from the BaBar model (squares). 



VI. CROSS CHECK FOR MODEL PREDICTIONS 

Using CP-tagged i^^^^Tr+vr^ samples, we can get q and c[ without any correlations, so 
the differences between q and provide a good test of our predictions on the differences 
discussed in Section II. The comparison between measured Acj = — Cj and the BaBar 
model predictions is shown in Fig. 9. We find good agreement between the data and the 
results obtained using the BaBar model, modified to account for the difference between 
K^sir+n- and KItt+ti-. 



VII. FINAL RESULTS AND IMPACT ON 7/^3 MEASUREMENT 

Our final results for Cj, Sj, c- and s'^ are shown in Table X and Table XI, respectively. 
The statistical uncertainties dominate for q and s,. The systematic uncertainty due to Aq 
and Asi which relate the strong phase difference of K^-k^-k^ and —>■ K^-k^-k" is 

comparable to — but does not dominate — all other contributions to the total systematic 
uncertainty. 

To see the impact of our results on the 7/03 measurement, we generate toy Monte Carlo 
D^K^ samples with 7/03 = 60°, 5b = 130° and = 0.1. The D^K^ sample 

is large enough so that the statistical uncertainty associated with B decays is negligible. 
We assume the reconstruction efficiency is 100% and that no background is present. We 
fit for 7/03, 5b, and 10,000 times by samphng Cj and Sj according to their uncertainties 
and correlations. We find the width of the resulting 7/^3 distribution, shown in Fig. 10, is 
about 1.7°. However, a small bias of 0.5° is observed, which is believed to be caused by the 
unphysical q and Si pairs (617 out of 8000) with + sf > 1. 

Comparing with a model uncertainty of 7° for BaBar [5] and 9° for Belle [6], great 
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TABLE X: Fit results for Cj and Sj. The first error is statistical, the second error is the systematic 
uncertainty (excluding Aq, Asj), the third error is the systematic uncertainty due to Acj and Asj 
that relate the Kgir'^ir^ and if£7r+7r~ Dalitz-plot models . 



Ci 



0.743 ± 0.037 ± 0.022 ± 0.013 0.014 ± 0.160 ± 0.077 ± 0.045 

1 0.611 ±0.071 ±0.037 ±0.009 0.014 ± 0.215 ± 0.055 ± 0.017 

2 0.059 ± 0.063 ± 0.031 ± 0.057 0.609 ± 0.190 ± 0.076 ± 0.037 

3 -0.495 ± 0.101 ± 0.052 ± 0.045 0.151 ± 0.217 ± 0.069 ± 0.048 

4 -0.911 ± 0.049 ± 0.032 ± 0.021 -0.050 ± 0.183 ± 0.045 ± 0.036 

5 -0.736 ± 0.066 ± 0.030 ± 0.018 -0.340 ± 0.187 ± 0.052 ± 0.047 

6 0.157 ± 0.074 ± 0.042 ± 0.051 -0.827 ± 0.185 ± 0.060 ± 0.036 

7 0.403 ± 0.046 ± 0.021 ± 0.002 -0.409 ± 0.158 ± 0.050 ± 0.002 



TABLE XI: Fit results for and s'^. The first error is statistical, the second error is the systematic 
uncertainty (excluding Ac,, Asj), the third error is the systematic uncertainty due to Acj and Asj 



i 


4 







0.840 ± 0.037 ± 0.023 ± 0.014 


-0.021 ± 0.160 ± 0.080 ± 0.036 


1 


0.779 ± 0.071 ± 0.039 ± 0.008 


-0.069 ± 0.215 ± 0.060 ± 0.047 


2 


0.250 ± 0.063 ± 0.029 ± 0.102 


0.587 ± 0.190 ± 0.072 ± 0.006 


3 


-0.349 ± 0.101 ± 0.057 ± 0.092 


0.275 ± 0.217 ± 0.067 ± 0.058 


4 


-0.793 ± 0.049 ± 0.029 ± 0.036 


-0.016 ± 0.183 ± 0.046 ± 0.042 


5 


-0.546 ± 0.066 ± 0.028 ± 0.038 


-0.388 ± 0.187 ± 0.056 ± 0.072 


6 


0.475 ± 0.074 ± 0.026 ± 0.081 


-0.725 ± 0.185 ± 0.065 ± 0.058 


7 


0.591 ± 0.046 ± 0.021 ± 0.011 


-0.374 ± 0.158 ± 0.059 ± 0.054 



improvement on the 7/^3 measurement can be achieved by using a model- independent ap- 
proach incorprating CLEO-c's results on the strong phase parameters q and Sj presented in 
this article. This will be realized at LHCb where using lO fb^^ of data a statistical error on 
7/03 of 5.5° is anticipated [7]. The weight of S ^ DK,D — > Kg7r~^TT~ in the combination 
of tree-level 7 measurements at LHCb, which is predicted to have sensitivity of 1° — 2° [21], 
depends upon the CLEO-c's results on the strong phase parameters q and Sj presented in 
this article. 

Sensitivity to New Physics is obtained through the comparison of 7/ 03 measured directly 
in tree-level processes and indirect determinations of 7/^3. One indirect determination, 
7/03 = (6714)°, arises from the intersection of the mixing and sin 2/3 contours in the (p, 
fj) plane [22]. The uncertainty is dominated by the LQCD calculations for mixing [23] and are 
expected to improve. Another determination of 7/^3 follows from the unitarity constraint 
7 = 180° — a — (3 = (7015)°. Here the uncertainty is dominated by the determination of 
a/0^ = (88;^)° from B 7nr,p7r,pp [22]. 
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FIG. 10: Toy Monte Carlo Fit Results for 7/(^)3. 

VIII. SUMMARY 

In summary, using 818 pb~^ of e^e" collisions produced at the V^(3770), we make a first 
determination of the strong phase pararneters, q and Sj, in Table X. From a toy Monte Carlo 
study with a large sample of B"^ D^K^ data generated with 7/03 = 60°, 5b = 130° 
and tb = 0.1, we find that the decay model uncertainty on 7/03 is reduced to about 
1.7° due to these new measurements. As a result, the precision of the 7/^3 measurement 
using — > D^K^ decays will not be limited by the strong interference effects in the 
Kg-K^n^ decay. The improved precision in the direct determination of 7/^3 enabled 
by this measurement of the strong phase parameters q and Sj enhances sensitivity to New 
Physics through the comparison with indirect determinations of 7/^3. 
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